Hubbard U and Hund's Exchange J in Transition Metal Oxides: Screening vs. 
Localization Trends from Constrained Random Phase Approximation 
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In this work, we address the question of calculating the local effective Coulomb interaction matrix 
in materials with strong electronic Coulomb interactions from first principles. To this purpose, 
we implement the constrained random phase approximation (cRPA) into a density functional code 
within the linearized augmented plane wave (LAPW) framework. 

We apply our approach to the 3d and 4d early transition metal oxides SrMOs (M=V, Cr, Mn) and 
(M=Nb, Mo, Tc) in their paramagnetic phases. For these systems, we explicitly assess the differences 
between two physically motivated low-energy Hamiltonians: The first is the three-orbital model 
comprising the t2 9 states only, that is often used for early transition metal oxides. The second choice 
is a model where both, metal d- and oxygen p-states are retained in the construction of Wannier 
functions, but the Hubbard interactions are applied to the d-states only ("d-dp Hamiltonian"). 
Interestingly, since - for a given compound - both U and J depend on the choice of the model, so 
do their trends within a family of these compounds. In the 3d perovskite series SrMOs the effective 
Coulomb interactions in the t2 g Hamiltonian decrease along the series, due to the more efficient 
screening. The inverse - generally expected - trend, increasing interactions with increasing atomic 
number, is however recovered within the more localized "d-dp Hamiltonian". Similar conclusions 
are established in the layered 4d perovskites series Sr2M04 (M=Mo, Tc, Ru, Rh). Compared to 
their isoelectronic and isostructural 3d analogues, the 4d perovskite oxides SrMOs (M=Nb, Mo, Tc) 
exhibit weaker screening effects. Interestingly, this leads to an effectively larger U on 4d than on 3d 
shells when a t^ g model is constructed. 

PACS numbers: 71.27.+a,71.10.Fd,71.15.Ap,71.45.Gm 
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I. INTRODUCTION 

Naively, the calculation of the Coulomb repulsion be- 
tween two charges may seem to be a simple textbook 
problem. However, calculating this repulsion for two 
charges in a solid is far from being trivial, as the elec- 
tronic polarizability screens the Coulomb potential lead- 
ing to a renormalized repulsion strength. Obtaining a 
quantitative estimation of the interactions, is particularly 
important in situations where the Coulomb interactions 
dominate the behavior of the system 1 . 

The description of the excitations of such strongly cor- 
related materials require a theoretical treatment that 
goes beyond the one-particle picture of band theory. Sim- 
ply identifying the Kohn-Sham spectra of density func- 
tional theory (DFT) in the local density approximation 
(LDA) 2 ~ 5 with the many-body spectra, for example, be- 
comes then an inappropriate oversimplification. Popu- 
lar methods for introducing many-body corrections into 
this description construct a multi-orbital Hubbard-type 
Hamiltonian, where explicit interaction terms of Hub- 
bard and Hund form are added. The resulting problem is 
then e.g. solved within a static mean field theory within 
the "LDA+U" scheme 6 or within dynamical mean field 
theory (DMFT) within the combined "LDA+DMFT" 
method 7-9 . However, the predictive power of such meth- 



ods crucially relies on a reliable assessment of the inter- 
actions. The ab initio calculation of these parameters is 
hence an important issue. 

Within the last years, tremendous progress has been 
made concerning this question, in particular with the 
advent of the so-called "constrained Random Phase Ap- 
proximation" (cRPA) as introduced in Ref. [10]. The 
cRPA provides a systematic first principles technique 
for the construction of low-energy Hamiltonians where 
not only the one-particle part of the Hamiltonian but 
also the interaction part is calculated from first princi- 
ples, that is without adjustable parameters. The start- 
ing point is the choice of an effective low-energy Hilbert 
space, and the cRPA aims at constructing the partially 
screened interaction that should be used as the bare inter- 
action within this space. The procedure can be viewed 
as the analogue for the interaction term of the famil- 
iar "downfolding" techniques used for the one-body part 
of the Hamiltonian 11,12 . Matrix elements of the effec- 
tive partially screened interaction are identified with the 
Hubbard and Hund's interaction matrices for the ef- 
fective low-energy Hamiltonian, and the corresponding 
parametrizations in terms of Hubbard U and Hund's rule 
J can be explicitly constructed. 

In this work, we have implemented the constrained 
random phase approximation (cRPA) 10 within the full- 
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potential augmented plane wave ((L)APW+lo) frame- 
work into the popular density functional code Wien2k 13 . 
Our approach allows for the calculation of the Coulomb 
interaction matrix elements in a localized basis set, and 
can be combined with many-body techniques working 
in this framework, such as the recent implementation of 
LDA+DMFT within the same basis set 14 . We illustrate 
the power of our approach on the 3d and Ad transition 
metal oxides SrM0 3 (M = V, Cr, Mn, Nb, Mo, Tc) and 
Sr 2 M0 4 (M = Mo, Tc, Ru, Rh). 



It has been argued early on 15 ' 16 that in systems where 
the behavior of the correlated orbitals is sufficiently close 
to a purely atomic description, the Hubbard and Hund 
interactions acquire the status of physical observables, di- 
rectly measurable e.g. in spectroscopy experiments. As a 
consequence of negligeable interorbital hybridizations, a 
natural basis set for the definition of U and J (the atomic 
one) emerges from the physical situation in these cases. 
Here, we adopt a more general point of view, defining 
Hubbard and Hund interactions even for materials where 
hybridization between states of different orbital charac- 
ter is too strong to allow for uniquely defined atomic 
orbitals to emerge. Even in this general case, Hubbard 
and Hund interactions can still be defined as useful aux- 
iliary quantities for many-body calculations, but are no 
longer direct physical observables. We illustrate this is- 
sue on the examples of 3d and 4d transition metal oxides 
(ranging from correlated metals to Mott insulators), for 
which we explicitly construct two different many-body 
models, demonstrating the dependence of U and J on 
this choice. We stress that, as a direct consequence of 
this issue, statements on values of U and J for a given 
material are not meaningful pieces of information unless 
supplemented by the complete description of the model 
to which they apply, both concerning the choice of the 
low-energy subspace and of the correlated orbitals. 



The paper is organized as follows: in section II, we 
discuss the context of our work with respect to previous 
methods for the determination of Hubbard and Hund in- 
teraction parameters. A general description of cRPA and 
of its implementation within the (L)APW+lo framework 
is given in section III, supplemented by technical issues in 
Appendix A and a discussion of a parametrization of the 
Hubbard interaction matrix in terms of Slater parame- 
ters in Appendix B. The results of cRPA calculations on 
the 3d and 4<i series of ternary transition metal oxides 
SrMOa with (M=V, Cr, Mn) and (M=Nb, Mo, Tc) and 
on the 4d layered perovskites Sr 2 M04 with (M=Mo, Tc, 
Ru, Rh) are described respectively in Sections IV and 
V. In Section VI, we discuss perspectives for fully first 
principles many-body calculations for correlated materi- 
als. We conclude with a brief summary of key results in 
Section VII. 



II. HUBBARD AND HUND INTERACTIONS: 
FROM SPECTROSCOPY TO FIRST 
PRINCIPLES CALCULATIONS 

A. Hubbard U : fitting parameter or calculated 
quantity? 

Historically, the Coulomb energy cost for placing two 
electrons on an atomic site within a solid - the Hub- 
bard interaction U - was introduced in the single-orbital 
Hubbard-Kanamori-Gutzwiller model 17-19 . U corre- 
sponds to the difference between the electron affinity and 
the ionization energy when respectively adding and re- 
moving an electron on the correlated shell (e.g. the d- 
shell) of a given atom: 

U = E(d n+1 ) + E(d n - 1 )-2E(d n ), (1) 

where E(d n ) is the total energy of a system for which 
n electrons fill a given d shell on a given atom. When 
the orbital character of the additional or missing charge 
is specified also, an interaction matrix U can be defined. 
Depending on the degree to which the correlated d-shell 
retains its atomic character even in the solid, its elements 
can be parametrized by a small number of (Slater) pa- 
rameters. 

An important step in the determination of U and J was 
the pioneering work of Sawatzky and co-workers 15,16 ' 20-22 
who fitted the multiplet structure measured in x- 
ray photoemission, absorption and Auger spectra with 
configuration-interaction cluster models involving a set 
of Slater integrals (F fe ) 23,24 . In a very simplified way, 
U = F can be obtained from the Mott gap whereas the 
other Slater integrals F k ,k > can be deduced by fit- 
ting atomic Hartree-Fock-like calculations to experimen- 
tal spectra. 

Attempts to calculate U from first principles typi- 
cally rely on constrained density functional approaches 25 . 
Such "constrained LDA" (cLDA) approaches 26-30 arc 
based on the observation that the energy of a system with 
increased or reduced particle number is in principle acces- 
sible within density functional theory. £/ cLDA is then de- 
fined as the derivative of the total energy while constrain- 
ing the occupation on a given shell. cLDA-type schemes 
have been implemented e.g. within the linear muffin-tin 
orbital (LMTO) framework assuming that the LMTO ba- 
sis functions used in a multi-orbital Hubbard model give 
a good representation of the localized degrees of free- 
dom. Later on, implementations of cLDA were worked 
out within various other electronic structure codes, for 
example, in the (L)APW+lo framework 31 , or in a basis of 
maximally localized Wannier functions (MLWF) 32 . We 
mention moreover the development of alternative linear 
response formalisms assessing the screened interaction by 
explicitly perturbing the system 33,34 . 

In 2004, Aryasetiawan and co-workers 10 proposed a 
scheme nowadays known under the name of "constrained 
Random Phase Approximation (cRPA)" . The cRPA is 
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an approximation to a systematic Wilson-like procedure 
to downfold a system to a low-energy Hamiltonian. 

The cRPA method has been implemented within the 
LMTO in the atomic sphere approximation (ASA) frame- 
work, employing the heads of the LMTOs as local or- 
bitals 10,35 and within the MLWF framework using the 
full-potential (FP) LMTO 36 and the (L)APW+lo ba- 
sis 37 . Materials- wise, the Coulomb interactions in several 
families of oxypnictides, parent compounds of the new 
Fe-based superconductors have been investigated within 
cRPA in the MLWF framework 38 . 3d, Ad and 5d transi- 
tion metals have been studied in Refs. [37,39,40]. Inter- 
estingly, even the pressure dependence of the Coulomb 
interactions becomes accessible within cRPA 41 . A recent 
attempt to go beyond the cRPA method uses an LDA+U 
electronic structure instead of the LDA one as input to 
the cRPA scheme 42 . This technique has been applied to 
NiO and elemental Gd. 

In recent years, several combined LDA+DMFT stud- 
ies have appeared, which use the ab initio Hubbard and 
Hund interactions calculated from cRPA. More specifi- 
cally, the spectral properties of the parent compounds of 
iron-based superconductors LaOFeAs 14 and FeSe 43 , and 
of the layered perovskite Sr 2 Ru04 44 and the spin-orbit 
Mott insulator Sr 2 Ir04 as well as the related compound 
Sr2Rh04 45 have been investigated. Except for the cases 
of Sr 2 Ir04and Sr 2 Rh04 these works however evaluated 
the Coulomb interactions in the MLWF basis set whereas 
the LDA+DMFT calculations were performed in a differ- 
ent basis set. 



B. Hubbard U: physical or auxiliary? 

As alluded to above, the definition of U as the dif- 
ference between electron addition and removal energies, 
directly accessible in spectroscopy measurements, sug- 
gested to consider the local Hubbard interaction as the 
central physical quantity characterising Coulomb inter- 
actions in the solid. A change of paradigm was triggered 
in 2003, when the proposal of a combined GW and dy- 
namical mean field scheme ("GW+DMFT") 48 ' 49 led to 
a generalisation of the notion of U, degradating it at the 
same time - at least from a conceptual point of view - to 
a purely auxiliary quantity. Indeed, in the GW+DMFT 
scheme, a dynamical interaction function Ui oca i(ui) is in- 
troduced and (in principle sclf-consistently) determined 
such that the resulting fully screened Coulomb interac- 
tion W takes it physical value. This is akin to repre- 
senting the density - within density functional theory 
- by an auxiliary (Kohn-Sham-) potential, or the local 
Green's function in dynamical mean field schemes by a 
local impurity bath Weiss function. Formally, this leads 
to the interpretation of screening as a two-step process: U 
is interpreted as the bare interaction within an effective 
subsystem with reduced number of degrees of freedom 
(in the case of GW+DMFT, a dynamical local impurity 
problem). The latter is such as to generate the physical 



screened Coulomb interaction W. In the same sense, the 
cRPA as proposed in Ref. [10], defines U by the require- 
ment that the corresponding fully screened interaction 



takes its physical value, and P sub is the polarization of 
the chosen subsystem. The original cRPA calculates both 
W and P sub from the random phase approximation, but 
the general treatment in Ref. [50] for example, makes 
it clear that this is not an essential ingredient. Other 
choices for P sub have been explored e.g. in Refs. [51, 52, 
53, 54]. 



C. Constrained Screening Approaches (CSA) 

In this paragraph, we give a unified description of 
the above type of approaches, which we will refer to as 
"constrained screening approaches" (CSA). They have 
in common to consider Hubbard interactions U as auxil- 
iary quantities, constructed in such a way that the phys- 
ical fully screened Coulomb interaction W of the origi- 
nal Coulomb problem in the continuum, possibly within 
some approximation, coincides with the fully screened 
Coulomb interaction of a chosen subspace when U is ap- 
plied as a bare interaction to this subspace. The first step 
is thus the choice of a subsystem with reduced number 
of degrees of freedom. In the "GW+DMFT" approach 
or in the approach of Ref. [54] which directly constructs 
an impurity problem, this is a local problem. In the con- 
strained random phase approximation (cRPA), as pro- 
posed in Ref. [10], this is a low-energy subspace for which 
a Hubbard model is to be defined. Other choices, as for 
example subspaces of reduced dimension are also possi- 
ble; see e.g. Ref. [52] for dimensional downfolding to a 
2d-lattice model or Refs. [53,54] for the direct construc- 
tion of a local model. 

The requirement of the physical W be represented by 
this effective problem (cither its local part only, as in 
GW+DMFT, or its full momentum-dependence as in the 
case of the lattice cRPA) leads to a formal relation of U 
and W involving the polarization P sub of the subsystem: 
U^ 1 - W^ 1 = P sufc ; that is, U is the partially screened 
Coulomb interaction, screened by all processes present in 
the original system except those contained in p sub . For 
the GW+DMFT scheme, one simply has P sub = P lm P, 
the polarization of the dynamical impurity model. In 
constrained screening approaches geared at the construc- 
tion of a lattice model (rather than a local one) one sets 
psub e q ua i to the polarization P d of a low-energy corre- 
lated subspace. For obvious reasons, in practical calcu- 
lations equivalent approximations should be made on W 
and P. The constrained random phase approximation 10 , 
can thus be viewed as a special case, in which both, W 
and P sub are calculated within the random phase approx- 
imation. 
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In practice, we start by choosing a set of degrees of 
freedom around the Fermi level in order to generate a 
correlated subspace C of the full Hilbert space. The con- 
ceptual division of the Hilbert space into a reference space 
C and the remaining degrees of freedom leads to the fol- 
lowing decomposition of the response of the system to an 
external perturbation, expressed by the polarization P: 



p _ psub _|_ p 



(3) 



Here, P sub denotes the polarization within the correlated 
subspace C, and, by definition, P r = P — P suh is a con- 
strained polarization, in which the contributions of the 
target degrees of freedom have been projected out. The 
constrained polarization leads to the partial dielectric 
function e r 

e r (l, 2) = 6(1 - 2) - J d3 P r (l, 3M3, 2). (4) 

We use here a shorthand notation that consists in rep- 
resenting time and possible spatial degrees of freedom 
by a number. In the case of the usual lattice cRPA, for 
example, this number would thus represent (r r) , in the 
case of local subsystems it would be the time variable r 
only. The partially screened interaction W r can then be 
defined as follows: 



W r (l,2) = J d3e r T 1 (l,3)w(3,2). 



(5) 



A simple algebraic manipulation shows that further 
screening of W r by the polarization p sub } or, equiva- 
lently, by the low-energy effective dielectric function 

e au6 (l,2) - 5(1-2)- J d3P sub (l,3)W r (3,2). (6) 

allows to retrieve the fully screened interaction W: 

W(l,2) = J d3e^(l,3)W(3,2). (7) 

This property of W r suggests that it can be identified 
as the effective bare interaction within the low-energy 
subspace. More precisely, matrix elements of the static 
value of W r in the localized basis set can be interpreted 
as forming the interaction matrices : 



U.... = (..\W r 



(8) 



The following remark is in order here: in fact, the ob- 
tained interactions are frequency- dependent because of 
the frequency dependence of the polarization P r (Eq. 5). 
This is a physical effect, since the response of the elec- 
trons in the solid to an external perturbation depends on 
the frequency of the latter. In particular, the electrons do 
not respond to a high-frequency oscillating electric field 
whose frequency exceeds any electronic energy scale, that 
is, the electronic polarizability vanishes at high frequency 
and the partially screened Coulomb interaction W r co- 
incides, in this limit, with the bare Coulomb interac- 
tion. cRPA in principle provides us with the spectrum of 



the dynamical interaction whose combination with many- 
body solvers such as extended dynamical mean field the- 
ory is currently receiving much attention 55-57 . In par- 
ticular, it was shown how inclusion of the high-energy 
tail of the energy-dependent Coulomb interaction leads, 
in the spectral functions, to pronounced shifts of weight 
to higher energies and a concomitant additional renor- 
malization at low energies 55,56 . Recently, an explicit ex- 
pression for these renormalisations was derived 57 , yield- 
ing a prescription for the construction of a low-energy 
Hubbard-like model with static Hubbard and Hund in- 
teractions where the high-energy screening effects are al- 
ready incorporated. Specifically, it was argued that the 
static interaction matrices to be used in such a construc- 
tion are precisely the ones obtained from constrained 
screening approaches in the static limit (w — > 0). These 
are the subject of the present work. 



III. 



THE CONSTRAINED RANDOM PHASE 
APPROXIMATION 

A. Constrained Polarization 



While the conceptual construction of the Hubbard in- 
teractions is general and in principle independent of the 
specific choice of the subsystem, the construction ac- 
quires its sense from a Wilson-like argument, when the 
reference system is a low-energy subspace: P sub = P d , 
where the subscript d refers to a low-energy correlated 
target space. The constrained random phase approxima- 
tion then corresponds to the choice of calculating both, 
the polarization of this subspace and the total polariza- 
tion within the RPA approximation. If we assume that 
the low-energy bands that span the correlated subspace 
do not energetically overlap with the remaining bands, 
the C-restricted polarization reads 

occ unocc 

P d (ry-u) = Y J E r dk (r)^(rW dl y(r')^ k (r') 

k,d k',d' 
1 



X 



1 



w - erf'fe' +^dk + irj oj + e d , k > - e dk - in 



(9) 



where e dk are the energies of the Bloch states, ip d k(r), 
that span C. The polarization P r is then deduced from 
Eq. (3) with P sub = P d . 



B. Definition of the correlated subspace 

The choice of the "target" degrees of freedom - or low- 
energy subspace - is obviously not unique. Two different 
"downfoldcd" Hamiltonians will not lead to the same ef- 
fective interactions. However, both choices should yield 
the same results for physical observables at the end, un- 
der the condition that both models are appropriate for 
catching the relevant physics of the system. 
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In the following, we assume that the band structure 
- as in the perovskite materials considered in this pa- 
per - is such that the orbital character of the "correlated 
species" (here the 3d or Ad orbitals) is dominantly spread 
over a subset of bands that are not entangled with bands 
of other orbital characters. In other words, the Hilbert 
space can be split into a correlated and an itinerant sub- 
space that - for every k-point in the first Brillouin zone - 
are separated in energy. The situation of entangled sub- 
spaces requires specific treatments that will be reported 
elsewhere (see also Refs. [37,40,53]). 

Our implementation is based on the full-potential 
(L)APW+lo code Wien2k 13 , and the construction of lo- 
calized orbitals according to the procedure implemented 
in Ref. [14]. In the case of separated bands in the sense 
described above, this construction results in a set of Wan- 
nier functions {l^^)} that span the "target" space C. 
The index m denotes an orbital quantum number, a is 
an atomic index inside a given unit cell centered at R 
and cr is the spin degree of freedom. 

The practical construction of the orbitals {l^^)} pro- 
ceeds as follows: one starts by choosing the dominant 
orbital character of the states that are to be considered 
as correlated, as well as an energy window W that con- 
tains at least all bands whose majority character is the 
correlated one. One then uses a set of localized atomic 
orbitals of the chosen character {|Xfcm)i and calculates 
their projection onto the Bloch states that lie in the en- 
ergy window W: 

\xt a m ) = £<^„lxSTM„>- (io) 

The set of wave functions {|Xfcm)} * s then subject to an 
orthonormalization procedure, yielding the orthonormal 
Wannier basis It is an easy task to verify that 

this set is indeed complete within the selected subspace. 

In order to be able to transform back and forth be- 
tween the Kohn-Sham and the Wannier basis, a set of 
"Wannier" projectors P^(k) arc defined as follows : 

P^Ak) = (<OV>L>, (ii) 

and with simplified notations, L = (m,a,a), 

\<t>HL) = ]T P£„(fc)l</W>- (12) 

The atom-centered Wannier functions in which the Hub- 
bard model is formulated are finally obtained from a sim- 
ple Fourier transformation from the reciprocal to the di- 
rect space: 

k 

where Af corresponds to the number of k vectors inside 
the first Brillouin zone. 

As mentioned before, in general the choice of the differ- 
ent subspaces is not unique. Depending on the material 



and the physical quantities considered, one has to decide 
for example whether to include ligand states in the model 
or not. For the oxide compounds considered in this work, 
two natural choices appear. We use the notations previ- 
ously introduced in Ref. [58]. 

"i2 S -^2g Hamiltonian" : Due to the octahedral crys- 
tal field of the oxygen ligands, the d states of the metal 
ion M + in SrM03 are split into t 2g and e g bands. The 
t 2g bands are partially filled (with one electron in SrV03 
three electrons in SrMnOs) whereas the e g states are 
empty. In a quite intuitive way, one can choose the 
t2 g bands as the low-energy degrees of freedom in the 
cRPA downfolding procedure. The correlated subspace 
C t2g thus includes only the t 2g degrees of freedom. To con- 
struct this subspace explicitly, a set of Wannier functions 
with the t2 g character is constructed out of Kohn-Sham 
bands within the energy window W t2g ■ The energy win- 
dow W t2s is chosen such that only t 2g bands are included. 
To calculate the constrained polarization P r , transitions 
from and to the target t 2g bands included in Wt 2s are 
removed from the total polarization (Eq. 3). A method 
based on the Kohn-Sham indices or on an energy win- 
dow can be employed to label the transitions that have 
to be eliminated, since the t 2g bands do not energetically 
overlap with p and e g ones. 

We finally obtain the interaction parameters that cor- 
respond to a C t2g -restricted lattice Hamiltonian. We call 
this model the u t 2g -t 2g Hamiltonian" , where the first in- 
dex refers to the majority character of the bands that are 
cut out from the screening process when going from P to 
P r , and the second refers to the majority character of the 
bands that are used to construct the Wannier functions 
in which matrix elements are taken. 

"d-dp Hamiltonian": Alternatively, one can choose 
as correlated subspace Cd p the space that contains also 
the e g degrees of freedom in addition to the t 2g ones. In 
an octahedral crystal field, however, the (dzz 2 -r 2 Ax 2 -y 2 ) 
orbitals strongly hybridize with the ligands, forming 
bonding and anti-bonding bands. The former are dom- 
inantly of ligand-p character, with substantial weight of 
the transition metal d contribution, however, while the 
latter are dominantly formed by the e g states, with ad- 
mixture from ligand p states. This situation suggests to 
use for the construction of the Wannier functions an ex- 
tended dp energy window Wd P that includes the whole d 
manifold as well as the p one. It could thus appear nat- 
ural in this context to construct a low-energy effective 
Hamiltonian, where both, d- and p-states, are treated 
as correlated states. Treating Coulomb interactions be- 
tween the in general rather extended p electrons as local, 
however, seems to be a more drastic approximation than 
treating these interactions in a static mean field fash- 
ion. In practice, it is thus more common to describe 
oxides within a hybrid model, where the Wannier func- 
tions are constructed within an extended dp window, 
while applying Hubbard interaction terms only to the 
rf-manifold. Analogous constructions have for example 
been performed for iron pnictide compounds 58 where the 
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Fe-d states were considered as correlated while "uncorre- 
lated" pnictogen or oxygen states had been included for 
the sake of the construction of sufficiently localized Wan- 
nier functions. The philosophy that is currently taken in 
this context is then to construct the constrained polar- 
ization P r from the total polarization by cutting off the 
transitions from and to the target d bands. The low- 
energy subspace described by the Hamiltonian is how- 
ever composed of the full d- and p-manifolds. Following 
the notation of reference 58 we call this hybrid model the 
ll d-dp Hamiltonian" . 



C. Calculated quantities 

1. Definition of matrix elements 

The cRPA method allows to calculate the effective in- 
teraction matrix elements of a low-energy lattice Hamil- 
tonian expanded in the localized basis set, {\4>rl)} ■ 



cubic harmonics, m m , includes in principle not 

only density-density interaction terms. Reduced inter- 
action matrices, U aa , U aa , for the density-density inter- 
action with cubic symmetry between opposite spin and 
parallel spin electrons respectively, can be defined as fol- 
lows : 



(14) 

= J J drdr'<l)* RlLl {r)4> R3L3 {r)W r {ry;uj) 

x0Wr>K 4 L 4 (r'). (15) 

By expanding the localized orbitals into the Kohn-Sham 
states within the energy window W, the interaction ma- 
trix elements read (Appendix A) : 

jjR 1 R 2 R 3 R4,r u \ _ j\f-2 e i(fei-Ri-fe 3 -R3+fe2-R2-fe4R4) 

iWfei)[iW*3)]* 



k±k2 k^k^ 



V\ V-2 V3 va £W 



(17) 



and the bare interaction, 

t RiR 2 R 3 R& 



V 



L\L 2 L 3 L^ 



R 1 L 1 <PR 2 L 2 \V\(pR 3 L 3 (PR i L l 



(18) 



While the procedure is in principle directly applicable 
also to non-local interactions, we will restrict the discus- 
sion in the following to local (that is intra-atomic) inter- 
actions (Ri = R 2 = -R3 = Ri). The interaction matrix 
elements are denoted 17f im!ra3m4 where m refers to an- 
gular quantum numbers indicating the dominant orbital 
character of the local basis set and S to the symmetry of 
the crystal field used for constructing the localized basis 
(in our case, cubic harmonics or spherical har- 
monics). The four- index interaction matrix within the 



-aa 
mm' 



cubic 
mm'mm' 
cubic 
mm'm'm 
JJ(J(J = Tycubic 
^ mm' — ^ mm'mm' 



= u: 
= u, 



= (^m<t>m'\W r (0)\<f> m ^ m ,} (19) 
= (Mm'\W r (0)]^,^) (20) 
— J mm' ■ (21) 



For practical reasons, and in order to simplify its fur- 
ther use within many-body techniques, it is convenient to 
parametrize these matrices by a reasonable number of pa- 
rameters. Common practice and physical insight suggests 
the popular choice of the Slater parametrization. Based 
on the assumption of spherical symmetry, this form is 
strictly valid only for atoms, but yields - for sufficiently 
localized atomic-like states - a good approximation for 
the correlated shells of a solid. 



2. Slater integrals 

The parametrization by Slater integrals 23 ' 24,59 makes 
use of the specific form of the angular momentum in- 
tegrals in the atomic case, reducing the problem to the 
calculation of a restricted number of radial integrals (3 
in the case of a d shell, 4 in the case of a / shell). Use- 
ful relations are summarized in Appendix B. The Slater 
integrals F fe (u;) can be employed for building the sym- 
metrized interaction matrix U^^ m mi {w) with cubic 



symmetry as follows: 



V, 



cubic 

raim^mzmA v 



x(tpk 1 u 1 tpk 2 ^\W r \u)\ll)k 3 v 3 i>k iV4 )[PL iVi {ki)]* PL 2 „ 2 (fc 2 )- 

(16) 

Analogously, one can compute the fully screened interac- 
tion, 



( w ) ^ ^ ^m 1 m' 1 S m2m ' 2 

m'^m'^m'^m'^ 

J2 a k ( m \,m> 2 , m' 3 ,m' 4 )F k ^))sJ 



fe=0 



(22) 



where a k are the Racah-Wigner numbers (Appendix B) 
and S is the transformation matrix from spherical har- 
monics to cubic ones. 

Different conventions exist for the Hubbard U param- 
eter. We follow the convention of identifying U with the 
Slater integral F° (Refs. [6,60]) related to matrix elements 
of the interaction as a simple average over all possible 
pairs of orbitals: 



U = 



j 21 

— ^ ttl (m,m',m,m')F' (23) 

(24) 



{21 + 1) 



mm' k—0 



Hund's exchange J is given by 



60 . 



f~Tij 2 £"feKm»)F*. (25) 



2Z(2Z + 1) 



m^m' k 
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Using Slater integrals manipulations, it can be shown 
that : 



')/14 



(F + 

(286F 2 + 195F 4 + 250F 6 )/6435 



(for d shells) 
(for / shells), 
(26) 

Analogously, one can define bare parameters, v and 
Jbare> and fully screened ones, W and J scr ccnod, using 
the appropriate Slater integrals. The underlying assump- 
tion of this construction is that all three, bare, partially 
screened and fully screened interactions are spherically 
symmetric, that is keep the form imposed by the purely 
atomic angular momentum integrals. In general, this can 
be expected to be a good approximation for the bare in- 
teraction if the orbitals are sufficiently localized and thus 
atomic-like. In how far it is a good description for the 
partially and fully screened interaction depends on in how 
far screening acts as to break the spherical symmetry, in- 
ducing additional anisotropies through hybridizations in 
the local environment in the solid. 

Based on Eq. 22, one can extract Slater- symmetrized 
effective interactions for the whole d shell, corresponding 
to density-density interactions within the t 2g subspace, 
the e g subspace and also between the t 2g and e 9 24 . For 
further comparison with the interactions between t 2g or- 
bitals but within the t 2g -t 2g model, we will restrict our- 
selves to the Slater-symmetrized interactions for the t 2g 
subspace only : 



Umm' 
Jm. 



f° + If 2 + — f 4 

49 49 

F° - A F 2 - — F 4 
49 441 



(27) 
(28) 



U„ 



,) = ±F 2 + — F 4 . (29) 
' 49 441 V ' 



3. Hubbard- Kanamori parameters 

In the case of large crystal field splittings, the assump- 
tion of spherical symmetry of the correlated states is no 
longer valid, and the appropriate parameters for con- 
structing the £2 S -restricted lattice Hamiltonian (t 2g -t 2g ) 
are the Hubbard-Kanamori parameters 24 . They are di- 
rectly calculated from the interaction matrix elements in 
the C t2g subspace as follows : 



the bare and fully screened interaction matrices, « cublc 
and W cublc , respectively. 



4- Calculated quantities: summary of notations 

The Hubbard interaction matrix elements with cubic 
symmetry are denoted U^^ 2Tn3m4 (Eq. 16), and the 
corresponding reduced interaction matrices, U 1 
(Eqs. 19 and 21). 

Within the d-dp Hamiltonian, the average interaction 
matrix elements between t 2g orbitals that are directly 
calculated with cubic symmetry are denoted U mm , U mm > 
and J m for the on-site intra^g, inter-i 2s and exchange 
interaction, respectively. Analogous quantities, but for 
the bare situation, v. 



aa JJ aa 

mm' ' mm' 



n, Jm C > and for the fully screened 
situation, W mm , J^ reened , can be introduced. 

We deduce from the Slater integrals parametrization, 
the Hubbard parameter U = F (Eq. 24) and Hund's ex- 
change J = (F 2 + F 4 )/14 for the d shell (Eq. 26). They 
are the parameters that should be considered for con- 
structing the interaction Hamiltonian downfolded into 
the d-dp low-energy subspace chosen (Eq. 22). Analo- 
gous parameters but within bare (unscreened) repulsions, 
v and Jbaro, or within fully screened repulsions, W and 
^screened, are introduced. 

Employing Eqs. 27, 28 and 29, Slater-average on-site, 
U mm , and exchange, J m , interactions between t 2g or- 
bitals within the d-dp Hamiltonian, can be extracted - 
and analogously for the bare, Vmm,!}^™ and the fully 
screened W mm , J^ rccncd cases. They correspond to a 
part only - the t 2g one - of the d-dp low-energy Hamilto- 
nian and therefore, should not be taken for parametrizing 
this Hamiltonian. 

Within the t 2g -t 2g Hamiltonian, the Hubbard- 
Kanamori terms, U,U',J (Eqs. 30, 31 and 32), re- 
fer to the interactions between t 2g orbitals within the 
f2 S -restricted Hamiltonian. U corresponds to the intra- 
orbital Coulomb repulsion, whereas U'{= U — 2 J with 
cubic symmetry) is the inter-orbital interaction which is 
reduced by Hund's exchange, J . They are the appropri- 
ate parameters for mapping the low-energy Hamiltonian 
downfolded into the t 2g subspace. Analogous parame- 
ters but within bare repulsions, V,j7barc, or within fully 
screened repulsions, W, J r scrC cnod, are also introduced. 



w 



1 



AT=3 



u = N E ^ 



cubic 
J\J Z w rnmmm 



m—1 



N=3 



rr 1 \ A rrcubic 

~~ N(N-l) U mm'm' 



(30) 



(31) 



m^m' 
N=3 



T _ I) E/ U mm r mm r, (32) 



N(N — 1) 



m^m' 



Analogously, one can define bare parameters, (V, J7bare), 
and fully screened ones, (W, j7" S creened), by considering 



IV. HUBBARD INTERACTIONS IN 
PEROVSKITE TRANSITION-METAL OXIDES 

SRM0 3 (M = V, CR, MN, NB, MO, TC) 

A. General trends 

We first consider the ternary transition metal oxides 
SrM0 3 (M=V, Cr, Mn, Nb, Mo, and Tc) in an undis- 
torted cubic perovskite structure. This is the stable 
structure for SrVC>3, while for the other compounds this 
is an idealisation. In this structure, the transition metal 
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Figure 1: (color online) Electronic band structures (top) and projected density of states (bottom) of SrMC>3 (M = V, Cr, Mn 
from left to right) obtained from LDA in the paramagnetic phase. The 3d t2 g states are highlighted in red (dash) , the 3d e 9 in 
blue (dash-dot) and the oxygen p states in maroon (solid). 



Table I: Lattice parameters used for cubic perovskites SrM03 
and energy windows W (in eV) for the d-dp and the ti a -ti a 
models, d and t2 g Wannier-like functions are constructed out 
of the Kohn-Sham states included in W. Because of the en- 
tanglement of the e g states with Sr-like d states, the d-dp 
model is not considered for Ad SrM03. 





a (A) 


w dp 


w t2o 


SrVOs 


3.842 


[-7.5,5.5] 


[-1.8,1.8] 


SrCrOs 


3.820 


[-7.5,4.7] 


[-1.7,1.0] 


SrMnOs 


3.805 


[-7.5,4.2] 


[-1.7,1.0] 


SrNbOa 


3.997 




[-3.0,2.8] 


SrMoOg 


3.976 




[-3.0,2.0] 


SrTcOa 


3.950 




[-2.6,1.3] 



ion is octahedrally coordinated with oxygen ligands, in- 
ducing a splitting of the d orbitals of the transition metal 
into tig and e g , where the e g form bonding and anti- 
bonding states with the oxygen p. The ti g states form 
a relatively narrow partially filled band, separated in en- 
ergy both from the e g and the p states. There are two im- 
portant factors that characterise trends in the electronic 
structure of transition metal (TM) perovskites: 1) the lo- 



calization of TM-d orbitals, which increases as the atomic 
number increases within the same period of the periodic 
table, and decreases as the atomic number increases in 
the same group, and 2) the hybridization between tran- 
sition metal d and 0-2p states, which depends on the 
localization of the d orbitals as well as on the energetic 
splitting between d and 0-2p, the charge transfer energy 
and ligand field A p( j. The interplay of these two factors, 
together with the filling of the d manifold, gives rise to 
extremely intriguing physics and chemistry. The DFT- 
LDA band structures of the 3d and Ad transition metal 
perovskites in their paramagnetic phases are shown in 
Figs. 1 and 2, respectively. The lattice parameters used 
in the calculations are listed in Table I. 

Due to the contraction of the d orbitals with increas- 
ing atomic number, the hybridization between d and p 
orbitals lessens from the early to the late oxides. Simi- 
larly, the charge transfer energy between d and p bands, 
A p d, also decreases, in qualitative agreement with optics 
experiments 61,62 . On the other hand, the charge transfer 
energy is larger in Ad than in 3c?, due to the larger orbital 
extension of the former (Fig. 2). 

The narrowing of the 3c? bands around the Fermi level 
makes the corresponding oxides prototypical for the in- 
terplay between electronic itineracy and localization 1,63 . 
The 3d-based materials are usually considered as more 
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"correlated" than their 4c? analogues due to the more lo- 
calized character of the 3c? orbitals. Within the Zaanen- 
Sawatzky- Allen classification 64 and the evolution of A p( j 
through the series, the early transition metal oxides are 
more prone to Mott localization whereas charge transfer 
physics is dominant in the late ones. As commonly found 
in the literature - although never calculated yet - the lo- 
cal intra-orbital interactions are believed to increase with 
the atomic number in a similar way as in atomic systems. 
We show below that the trend for these interactions can 
differ from the one for atoms and we rationalize this dis- 
crepancy within the antagonism between the effects of 
the orbital localization and the screening. 



B. 3d perovskites: experimental facts 

We now give a brief - and by no means exhaustive - 
summary of experimental results on the oxide pcrovskitc 
families, SrM0 3 (M= V, Cr, Mn). 

SrVOs: Among the chosen oxides, SrV03 has been 
intensively investigated both experimentally and theoret- 
ically as a prototypical TM oxide that falls in the inter- 
mediately correlated regime. Experimentally it is found 
to be an undistorted paramagnetic metal which has been 
well characterized by experiments like optics, thermody- 
namical measurements, transport or angle-integrated and 
angle- resolved photoemission spectroscopy 1,65-69 . Pho- 
toemission spectra (Refs[ 67,68]) show a lower Hubbard 
band at about 1.7 eV binding energy and a quasipar- 
ticle peak characterized by the renormalization factor 
Z k> 0.6. An electron addition peak d 1 — > c? 2 has been 
identified by inverse photoemission at about 2.5 — 3 eV 69 . 
DFT-LDA can give a qualitatively correct description of 
the band structure of SrVC>3, as shown in Fig. 1, but a 
quantitative description including the band renormalisa- 
tion requires a more accurate treatment of the correlation 
effects in the c?-manifold 14 < 52 ' 65 < 70 . 

SrCr0 3 : Because of difficulties in the synthesis of 
SrCrC>3, only few studies have been conducted. Early 
work on single crystals with cubic structure 71 shows a 
metallic behavior with a Pauli paramagnetic susceptibil- 
ity. This is in disagreement with a more recent study 
on polycrystals reporting transport, thermal conductiv- 
ity and magnetic susceptibility in favor of a non-magnetic 
insulating state 72 . SrCrOs crystals have been also re- 
cently reinvestigated within x-ray diffraction studies 73 . 
The common belief found in the literature supported 
by density functional calculations - is that a structural 
transition from a non-magnetic orbitally-degenerate cu- 
bic to a distorted tetragonal- maybe antiferromagnetic- 
structure with orbital ordering could appear below 70 

j^73-75 

SrMnOs: At room temperature, SrMnOs is found 
in an hexagonal phase 76-79 but a cubic phase can be 
quenched and stabilized in a metastable state down to 
low temperatures. Both phases are deeply insulating. A 
G-type antifcrromagnetic ordering with a magnetic mo- 



Table II: (top) Slater integrals (in eV) for the d-dp model of 
SrMOa (M=V, Cr, Mn) corresponding to screened (W r ) and 
bare (v) Coulomb interaction, (bottom) Slater-symmetrized 



effective interactions (in eV) between t2 g orbitals. 


(eV) F° F 4 /F 2 J 


Fbarc F /F |barc ^bare 


SrVOs 3.2 0.795 0.85 
SrCrOs 2.9 0.781 0.85 
SrMnOs 2.8 0.774 0.89 


19.5 0.652 1.06 

20.1 0.628 1.06 

21.2 0.625 1.11 


Umm ^ mm' J m 


v v , 7 bare 


SrV0 3 4.1 2.8 0.65 
SrCrOs 3.9 2.6 0.65 
SrMnOs 3.8 2.4 0.68 


20.7 19.1 0.81 

21.4 19.7 0.82 

22.5 20.8 0.86 



ment around 2.6 [Ib emerges below the Neel temperature 
T N ~ 260K 77 . According to x-ray photoemission and 
absorption 80-82 , the spectroscopic properties of SrMnOs 
are rather involved. Indeed, the O-p states are strongly 
entangled with the lower Hubbard band of t 2g character, 
and the conduction band has been proposed to be of e g 
character 80 . These facts strongly question the validity 
of a pure t 2g model for the description of the low-energy 
spectra. We nevertheless present both, a t 2g and d-dp 
model, for the sake of comparison with work done in the 
literature and assessing trends along the perovskite se- 
ries. 



C. Hubbard parameters within the d-dp 
Hamiltonian 

We first consider the d-dp low-energy Hamiltonian for 
the 3c? series. The Slater integrals (Appendix B, Eq. 35) 
for the screened interactions, F fc , and for the bare inter- 
actions, F^ are are collected in Table II. The monopole 
part, F , is more efficiently screened than the multipolc 
parts, (F 2 ,F 4 ) 15 > 22 < 83 . The bare ratio F 4 /F 2 | barc is close 
to the atomic value around 0.63. In contrast, the par- 
tially screened ratio F 4 /F 2 deviates from this limit. This 
illustrates the importance of calculating the whole set 
of three Slater integrals for accurately parametrizing the 
four-index Hubbard interaction matrix, rather than de- 
ducing F°, (F 2 + F 4 )/14 from two independent relations 
and then assuming F 4 /F 2 = 0.63 as a third one. 

The Slater integrals can still be used in order to deduce 
effective interactions between t 2g orbitals (from Eq. 27 
to Eq. 29). The values of these interactions are shown 
in Table II. The validity of the Slater parametrization 
can be assessed by comparing these values to the ones 
obtained from a direct calculation with cubic symmetry 
(Appendix C). For SrVC>3, the direct calculation gives for 
the intra-t2g interaction U mm = 4.0 eV and J m = 0.60 
eV, hence in reasonable agreement with the values calcu- 
lated from the Slater integrals in Tab. II. For the SrM0 3 
series below, we will hence refer either to U mm , J m or 
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Figure 2: (color online) Electronic band structures (top) and projected density of states (bottom) of SrMC>3 (M = Nb, Mo, Tc from 
left to right) obtained from LDA in the paramagnetic phase. The id t2 g states are highlighted in red (dash), the id e g states in blue 
(dash-dot) and the oxygen p states in maroon (solid). 



In Fig. 3, the t^ g interactions for the bare, partially 
and fully screened cases, are shown. Along the 3d se- 
ries, the bare interactions v mm and J„ aro increase with 
the number of electrons in the d manifold. This is due to 
the increasing localization of the Wannier d orbitals from 
SrVC>3 to SrMnC>3 within the downfolded d-dp Hamilto- 
nian. A similar behavior is expected for hydrogenoid 
systems within the Slater rules, for which the atomic d 
radial extension decreases from V 4+ to Mn 4+ . The less 
extended the orbitals, the higher the Coulomb repulsion. 

The effective interactions, U mm and W mm , differs from 
the bare interaction v mm because of the screening that 
arises from the electronic polarizability and lowers the 
Coulomb repulsion. Interestingly, W mm decreases with 
the atomic number as a consequence of increasing screen- 
ing. Screening therefore counteracts the trend of increas- 
ing orbital localization that tends to increase the inter- 
action matrix elements towards the end of the 3d se- 
ries. The ratio W mm /v mm quantifies the increase of the 
screening from SrVC-3 to SrMnC>3. The smaller this ra- 
tio the stronger is the screening. We obtain 6.7/100 in 
SrV0 3 , 4.6/100 in SrCr0 3 and 4/100 in SrMn0 3 . 

Electronic screening is mediated by the creation of 
particle-hole and plasma excitations. At the RPA level, 
its strength varies with the inverse of the energy differ- 
ence between occupied and empty states. The knowledge 
of the DFT-LDA band structure (Fig. 1) thus allows us 



to analyse the efficiency of the different screening chan- 
nels. The charge transfer energy between the oxygen p 
and the transition metal d states decreases from SrV03 
to SrMn03. The Kohn-Sham states with ligand p char- 
acter go up toward the ti g ones. Furthermore, the crystal 
field splitting to the e g states decreases. While the abso- 
lute magnitude of the energy difference between p- and 
e g - states may depend on the approximation used in the 
band calculation (and in particular, is slightly underesti- 
mated in the LDA), the overall trend is not an artifact 
of DFT-LDA 1,62 . Experimentally, it has been evidenced 
e.g. in optical measurements 61 . As a consequence, dd 
and pd screening channels give stronger contributions in 
SrMnOs than in SrVC>3. However, these screening chan- 
nels are only part of the multiple screening channels that 
contribute to the reduction of the Coulomb repulsion 
from v mm to W mm . For example in SrVC>3, v mm = 20.7 
eV whereas W mm = 1.4 eV. As already pointed out in 
Ref. [39], the effects of the screening channels are not 
additive. Nevertheless, their contribution to the total 
polarizability can be estimated by evaluating partially 
screened interactions. 

The evaluation of the quantity U mm /v mm and the 
comparison with W mm /v mmi quantify the weight of 
the dd transitions in the total polarization. For ex- 
ample in SrVC>3, U mm /v mm equals 19.8/100 whereas 
equals 6.7/100. In SrCr0 3 , U mm /v mm = 
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(eV) (eV) 




I i i i I Lj i i_ 

SrV0 3 SrCr0 3 SrMn0 3 SrVO, SrCrO, SrMnO, 




v m m- Low-energy screening channels, such as dd or pd, 
are thus responsible for about half of the total screen- 
ing. On the other side, the transitions at higher energies 
also participate to the total polarizability leading to the 
strong reduction of the bare Coulomb repulsion. 

Hund's exchange J is deduced from the Slater integrals 
F 2 and F 4 (Tab. II). In the 3d SrM0 3 series, J goes 
from about 0.8 eV to Jbaro ~ 1.0 eV. The corresponding 
exchange interaction between t 2g , J m ~ 0.6 — 0.7 eV, is 
smaller than J and the values agree with the ones usually 
employed in models for such oxides. 

The fully screened exchange interaction J^ reened , cal- 
culated from the total polarization, can be evaluated for 
the series. J^ rccnod slightly decreases with the d electron 
number in contrast to the unscreened Hund's exchange 
J^ l are . Cutting off transitions from and to the d Kohn- 
Sham bands within Wt 2 , J m still increases from SrVO 3 
to SrMn03. In contrast to U mm , J m hence reproduces 
the atomic-like trend. The weaker dependence of the ex- 
change interactions on the screening results in J m varying 
in a similar way as J^ arc : . 



D. Hubbard parameters within the ti g -ti g 
Hamiltonian 



Figure 3: (color online) Middle panel: On-site Hubbard interac- 
tion U m m (left) and exchange interaction J m m (right) between t2 g 
orbitals within the d-dp (black curve) low-energy Hamiltonians for 
3q! SrM03, compared to the on-site U and J within ti g -ti g (red 
curve). In dashed line with open circles, the Hubbard parameter 
U = F° and Hund's exchange J = (F 2 + F 4 )/14 are shown. Top 
panel: Same but for the bare interactions between t2 g orbitals. Bot- 
tom panel: Same but for the fully screened interactions between t2 g 
orbitals. 



18.2/100 and W mm /v mm = 4.6/100, while in SrMn0 3 , 
U mm /vm m = 16.9/100 and W mm /v mm = 4.0/100. The 
partial screening without the dd channels, increases more 
slowly through the early 3d TM oxides than the to- 
tal screening. The interactions U mm thus correspond 
to an intermediate situation between v mm and W mm , 
where the increase in partial screening counterbalances 
the increasing localization of the orbitals. As a result, 
even though W mm decreases and v mm increases with the 
atomic number, U mm depends barely on it. 

The values of U m m indicate how strong the contribu- 
tion of the dd screening channels is to the total polariz- 
ability. These channels are only partly responsible for the 
screening that reduces v mm - The role of the pd screen- 
ing channels can also be estimated. By computing the 
partially screened interaction after the elimination of all 
pd transitions (as well as dd), one gets 11 eV in SrVOa, 
10.8 eV for SrCr0 3 and 10.7 eV for SrMn0 3 . Screening 
without these channels has thus reduced the interaction 
by roughly a factor two, compared to the bare interaction 



Within the ti g -ti g Hamiltonian, the ^-projected local 
orbitals within the energy window Wt 2g (Tab. I) lead to 
"extended" ti g Wannier orbitals. The charge transfer en- 
ergy and the hybridization between the ti g and the oxy- 
gen ligand p bands, are responsible for the hnite weight 
of the ti g Wannier functions on the oxygen atomic sites. 
The tail and hence the extension of the so constructed lo- 
cal orbitals, increases when the pd charge transfer energy 
becomes smaller, as happens for the 3d SrMOs series. 

The tig local orbitals are hence more extended in 
SrMn03 than in SrVOs, in contrast to the atomic or- 
bitals. It implies that i) the unscreened interaction 
V within the tig-tig model does not increase with the 
atomic number as v mm in the previous d-dp model (Fig. 3 
and Tab. II) and ii) the values of V are smaller than the 
values of v mm within d-dp. 

The interaction values that would be appropriate for 
the t2 9 -restricted Hubbard Hamiltonian are given in 
Tab. III. Their evolution along the series is illustrated in 
Fig. 3. A comparison with the cRPA values from the lit- 
erature but with different frameworks for constructing ti g 
local orbitals, is in order here: In SrV03, Aryasetiawan 
and co-workers 39 report U — 3.5 eV within the head of 
LMTO's, while Miyake and Aryasetiawan 36 give U = 3.0 
eV,J = 0.43 eV within MLWF. We obtain U = 3.2 eV 
and J = 0.46 eV. 

We now turn to the values through the early series of 
the 3d TM oxides. U significantly lessens from S1WO3 
to SrMn03 and the decrease of U is more pronounced 
than U mm within the d-dp model. The decrease of 11 
is still due to the screening which gets larger with the 
atomic number. More precisely, the p— > ti g and ti g —¥ 
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Table III: Hubbard-Kanamori U, bare V and fully screened W 
interactions (in eV) and corresponding exchange interactions, J , 
Jbaro an d ^screened between t^g orbitals within the t2 g -ti g down- 
folded Hamiltonian for the early 3d series SrMC>3 (M= V, Cr, Mn) 
and Ad (M= Nb, Mo, Tc). The inter-orbital interactions W coin- 
cide with U - 2 J. 



(eV) 


V 


Cr 


Mn 


Nb 


Mo 


Tc 


U 


3.2 


2.7 


1.8 


3.0 


3.0 


2.9 


J 


0.44 


0.42 


0.39 


0.29 


0.31 


0.31 


V 


16.1 


16.4 


16.2 


10.7 


11.6 


11.8 


jTbare 


0.55 


0.55 


0.53 


0.38 


0.40 


0.39 


W 


0.9 


0.4 


0.3 


0.9 


0.5 


0.4 


^/screened 


0.30 


0.17 


0.12 


0.24 


0.19 


0.16 


U/V x 100 


19.8 


16.4 


11.1 


28 


25.8 


24.6 


W/V x 100 


5.6 


2.4 


1.8 


8.4 


4.3 


3.4 



e g channels contribute more and more to the screening 
from SrV0 3 to SrMnOa, because of the p and e g Kohn- 
Sham bands that go closer and closer to the Fermi level. 
Quantitatively, the ratio U /V is about twice larger in 
SrV0 3 than in SrMn0 3 (Tab. III). The fully screened 
interactions, W, on the other hand, decrease through the 
series. 

In contrast to J m and J^ arc within the d-dp model, 
the exchange interactions, J and Jbaro within the t 2g -t 2g 
model slightly decrease with the number of d electrons. 
Considering the total polarization leads to a significant 
decrease of Screened through the series. Since the re- 
duction is weaker for J than for ^screened, the screening 
induced by the t 2g t 2g transitions must be responsible for 
the behavior of ^screened- As U, the exchange interactions 
therefore depend on the extension of the local orbitals: J 
is smaller within the t 2g -t 2g Hamiltonian than J m within 
d-dp. 



E. 4d perovskites: experimental facts 

We now turn to the 4d scries SrM0 3 (M= Nb, Mo, 
Tc). The lattice parameters used in the calculations are 
summarized in Table I. 

SrNb0 3 : This compound is usually obtained with a 
non-stoechiometric perovskite structure . When doping 
into Sr^NbOs with x > 0.80, a cubic perovskite phase is 
observed exhibiting a poor paramagnetic metallic behav- 
ior at temperatures below 300 K 85 . 

SrMo0 3 : It is, in contrast, an excellent paramagnetic 
metal- even the Ad transition metal oxide that displays 
the highest electrical conductivity 86 . With two electrons 
on the t2 g shell, it is an electronic analogue of SrRu0 3 
that has two holes on the t 2g s but larger correlation ef- 
fects due to the Van Hove singularity found in its density 
of states 44 . 

SrTc0 3 : Because of the radioactivity of technetium, 



this compound has been less studied, and only struc- 
tural and magnetic properties are known. A huge Ncel 
temperature of 1023 K - accompanied by a G-type anti- 
ferromagnetic ordering with the magnetic moment 2.1/ib 
below TV - has been recently discovered 87 . In particular, 
this high Tn is larger than in the 3d analogue SrMn0 3 . 
The large T/v in SrTc0 3 was first interpreted within den- 
sity functional calculations 87-89 , e.g. in terms of the 
larger covalency of the Tc-0 bonding compared to Mn- 
O 87 . Recently, another scenario has been put forward for 
interpreting the difference of magnitude in Tn between 
SrTc0 3 and SrMn0 3 , based on the proximity of SrTc0 3 
to the Mott transition in the presence of large Hund's 
exchange at half- filling 90 . 



F. Hubbard interactions in 4d perovskites 

In the Ad SrM0 3 series, the evolution of the Kohn- 
Sham bands from Nb to Tc oxides is analogous to the one 
in the isoelectronic and isostructural 3d compounds, but 
the trends are less pronounced (Fig. 2). The pd charge 
transfer energy decreases more slowly in the Ad series, 
and so do the crystal and ligand fields that split the t 2g 
and e g bands. As a result, the transition metal e g states 
overlap with strontium 5s /Ad. For this reason, only the 
t2 g Hamiltonian is considered here. 

The atomic-like character of the t 2g Wannier basis is 
evidenced by the atomic-like behavior of the bare inter- 
action V in the Ad series (Fig. 4): V becomes larger with 
the atomic number, as in 3d when considering the d-dp 
Hamiltonian. The local basis designed for the t2 g degrees 
of freedom, is atomic-like for 4c? oxides in contrast to the 
one for t 2g in 3d. The difference in Ad comes from the 
larger charge transfer energy which leads to a smaller 
tail of the downfolded orbital with t 2g character on the 
oxygen sites. 

However, the orbital extension of the t 2g local orbitals 
has to be larger in Ad than in 3d oxides since the bare 
interaction V are twice smaller in the former (Tab. III). 
This makes sense in an atomic-like basis in which the 
extension of the Ad atomic wavefunctions is larger than 
the extension of 3d. 

Since the pd charge transfer energy as well as the t 2g - 
e g splitting decrease slower in Ad oxides, the pd and the 
t 2g e g channels arc not as efficient in screening as in the 
3d analogues. This is quantitatively highlighted by the 
ratios U /W (Tab III), which are larger in Ad and decrease 
slower through the Ad series. 

Consequently, U in Ad is almost constant with the d 
electron number, like U mm in the d-dp model for 3d TM 
oxides, whereas U in 3d significantly lessens until becom- 
ing smaller in SrMn0 3 than in SrTc0 3 (Fig. 4). This is 
an effect of the screening which has stronger impact in 
the 3d extended Wannier basis than in the atomic-like Ad 
one. 

The exchange interactions in Ad oxides behave in a 
similar way as in 3d (Tab. III). The exchange interactions 



13 




♦ 



Hexagonal SrMnO^ 




1 1 


i 








4d 




1 



(1.3 - 
0.25- 
2- 
0.15- 



Figure 4: (color online) Middle panel: On-site Hubbard-Kanamori 
interaction tt (left) and exchange interaction J (right) between 
t2 g orbitals within the tig-t^g low-energy Hamiltonian for 3d (red 
curve) and 4d (indigo curve) SrM03 (see also Tab. III). Top panel: 
Same but for the bare interaction V and J/bare between tig orbitals. 
Bottom panel: Same but for the fully screened interaction W and 
i7screened between tig orbitals. 




Figure 5: (color online) DFT-LDA band structure of the four- 
layers hexagonal unit cell of SrMn03 in the paramagnetic phase. 
Lattice parameters a = 5.461A°,c = 9.093A° are taken from 76 . 



Table IV: Screened and bare Slater integrals (in eV) for the d-dp 
model in cubic and hexagonal SrMn03. 



(eV) 


F° 


F 4 /F 2 


J 


F? 

^ bare 


F /F [bare 


"/bare 


cubic 


2.8 


0.774 


0.9 


21.2 


0.625 


l.l 


hexagonal 


3.1 


0.777 


0.9 


21.9 


0.621 


l.l 



than in the cubic one. This is an effect of the Wannier or- 
bital localization rather than a screening effect since also 
the bare interaction v is enhanced. Hund's exchange J 
does not change much between the two crystal structures. 



in Ad are about 0.1 eV smaller than in 3d. This is due to 
the higher extension of the Ad Wannier orbitals which is 
also responsible for the smaller bare interactions V. 



V. HUBBARD INTERACTIONS IN THE 
LAYERED PEROVSKITE OXIDES SR 2 M0 4 (M = 
MO, TC, RU, RH) 

A. General features 



G. The influences of structural distortions: 
example of SrMnOj 



the 



At room temperature, SrMnC>3 crystallises in the four- 
layers hexagonal structure (P63/mmc) with a = 5.461 A, 
c = 9.093 A 76 . The DFT-LDA band structure is shown 
in Fig. 5. 

The four-times larger unit cell leads to a backfolding of 
bands in the first Brillouin zone. A set of i2g-like bands 
(Fig. 5) can be identified with the twelve bands around 
the Fermi level and a set of e s -like bands with the eight 
ones above. A d-dp model is constructed rather than a 
tig one. This allows for a direct comparison of the Slater 
integrals (Eq. 35) calculated within the d-dp model built 
for the cubic crystal structure (Tab. IV). 

U as well as v are slightly bigger in the hexagonal phase 



The second class of oxides that we have investigated 
are the layered perovskites Sr 2 M04 (M= Mo, Tc, Ru, 
Rh) where M 4+ is a Ad transition metal. The lattice 
parameters used in the electronic structure calculations 
for the paramagnetic phase are given in Table V. 

Sr2Mo04: This compound exhibits a metallic behav- 
ior over a wide range of temperature between 80 mK and 
300 K, with a resistivity increasing between 2 and 10 
mf2.cm . It is usually investigated under the possibility 
that it could exhibit analogous electronic properties than 
Si'2Ru04, although a superconducting state has not been 
reported down to 25 mK. 

Sr2Tc04: Due to the radioactivity of technetium ele- 
ments, only structural properties are known for Sr2Tc04. 
An undistorted layered perovskite structure is considered 
below. 
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Table V: Lattice parameters used for the layered Sr2M0 4 per- 
ovskitcs and energy windows W (in eV) for the d-dp and t2 9 -t2g 
models. 





a (A) 


c(A) 


w dp 


w t2o 


Sr 2 Mo0 4 


3.917 


12.859 


[-10,8.5] 


[-2.0,2.0] 


Sr 2 Tc0 4 


3.902 


12.720 


[-10,7.5] 


[-2.6, 1.3] 


Sr2Ru04 


3.863 


12.724 


[-10,6.9] 


[-3.0,1.0] 


Sr 2 Rh0 4 


3.854 


12.880 


[-10,6.0] 


[-3.3,0.5] 



Sr 2 Ru0 4 : The resistivity of this compound obeys a 
T 2 law below 30 K, evidencing a Fermi liquid behavior. 
A strong anisotropy of the transport properties - due 
to its layered structure - was reported 92 . A large mass 
enhancement and a low coherence scale have been deter- 
mined experimentally. Sr2Ru04 also exhibits an uncon- 
ventional superconductivity below 2 K 93 . 
The three-sheet Fermi surface determined experimentally 
is reasonably well described by DFT-LDA calculations, 
but the enhancement and the anisotropy in the mass are 
missed. The largest mass enhancement surprisingly ap- 
pears for the widest d xy band as determined by quantum 
oscillations experiments 93 . Several LDA+DMFT calcu- 
lations have been carried out to reproduce the low-energy 
properties of the photoemission spectra 44 ' 94-96 . 

Sr 2 Rh04: The symmetry is lowered from the K2N1F4 
class by 11° rotation around the c-axis of the RhOg 
octahedra 97 . It is a paramagnetic metal down to 
36 mK 98 . Spin-orbit coupling (SOC) was found to 
be relevant in addition to electronic correlations 99-101 . 
LDA+DMFT calculations with SOC have been recently 
performed for the distorted structure and compared to 
the isoelectronic but Mott insulating Sr 2 Ir04 45 . In par- 
ticular, it was shown how the interplay of SOC, corre- 
lations and structural distortions leads to a suppress of 
the effective orbital degeneracy, leaving only two orbitals 
at the Fermi level. In the following, for computational 
reasons, the undistorted crystal structure of Sr 2 Rh04 is 
considered. 

The DFT-LDA band structures for the paramagnetic 
phases are shown in Fig. 6. Bands with (d xy ,d xz ,d yz ) or- 
bital character emerge around the Fermi level, whereas 
the {d^ z 2_ r 2,d x 2_ y i) orbital character is found above and 
the oxygen p's lie below. The band with d xy character 
leads to a quasi-two dimensional Fermi surface whereas 
the degenerate bands with (d xz ,d yz ) character give rise to 
a quasi-one dimensional Fermi surface. In Sr 2 Mo04 and 
Sr 2 Tc04, the e g states are entangled with the strontium- 
like d's but come closer to the i 2 g's when the Ad elec- 
tron number increases. On the other hand, the pd charge 
transfer energy decreases with this number, since the p's 
go up toward the Fermi level. This decrease makes larger 
the screening from the p channels, in an analogous way 
to in the early transition metal oxide series. 






Figure 6: (color online) DFT-LDA paramagnetic band structures 
for layered perovskites Sr2M0 4 . (top): From left to right: M = 
Mo, Tc. (bottom): From left to right: M = Ru, Rh. 



B. Hubbard parameters within the d- 
Hamiltonian 



As previously, also for the early TM oxides within the 
layered perovskite structure Sr 2 M04, we consider first 
the d-dp Hamiltonian. The Ad Wannier orbitals are con- 
structed out of the Kohn-Sham bands within the energy 
window W dp (Tab. V). 

The Slater integrals are given in Tab. VI. The value 
of the ratio F 4 /F 2 |b a re is close to the one calculated 
for Ad atoms 102 . A significant deviation is obtained for 
the screened ratio of the Slater integrals, which is even 
stronger than the one in the 3d SrM03 series. This seems 
natural, given the anisotropy of the structure and the 
screening which increases with the orbital extension. 

A cubic approximation can be used for deducing a set 
of interactions between (d xy ,d xz ,d yz ) local orbitals from 
the Slater integrals (Eqs. 27, 28 and 29). The values are 
shown in Tab. VI and can be compared to the matrix 
elements calculated directly (Tab. VII and see also Ap- 
pendix D). The latter are anisotropic within the plane of 
the TM and oxygen octahedra. The Slater parametriza- 
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Table VI: (top) Static and bare Slater integrals (in eV) for the d-dp 
Hamiltonian in Ad Sr2M04. (bottom) Slater-symmetrized effective 
interactions between t2 g orbitals. 



(eV) 


F° 


F 4 /F 2 


J 






* bare 


F 4 /F 2 |barc 




Mo 


3.26 


0.862 


0.67 






14.50 


0.684 


0.86 


Tc 


3.19 


0.850 


0.70 






15.25 


0.673 


0.90 


Ru 


3.23 


0.838 


0.74 






15.97 


0.669 


0.94 


Rh 


3.44 


0.820 


0.78 






16.77 


0.663 


0.98 




Urnm 






.7 


n 






rbarc 


Mo 


4.0 


3.0 




0.50 


15.5 


14.1 


0.66 


Tc 


4.0 


2.9 




0.53 


16.3 


14.9 


0.69 


Ru 


4.1 


2.9 




0.56 


17.0 


15.6 


0.72 


Rh 


4.3 


3.1 




0.59 


17.9 


16.4 


0.75 
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Fi gure 7: (color online) Bottom panel: On-site interaction Umm 
(left) and exchange interaction J m (right) between t2 g orbitals 
within the d-dp (black curve) model for Sr2M04 and on-site in- 
teractions U (left) and {J (right) but within t2 g -t2g (red curve). In 
dashed line with open circles, the Hubbard parameter U = F° and 
Hund's exchange J = (F 2 + F 4 )/14 are shown. Top panel: Same 
but for the bare interactions between the t2 g orbitals. 



tion is more accurate for the late materials, Sr2Tc04 and 
Si'2Rli04, for which the spherical approximation of the 
Ad orbital is better due to the smaller hybridization with 
the ligands. 

The bare on-site and exchange interactions increase 
with the number of the 4c? electrons (Fig. 7), suggest- 
ing a rather atomic-like behavior of the Ad Wannier or- 
bitals within the d-dp Hamiltonian. This agrees with 
the atomic expectations based on the Slater rules, ie a 
higher localization due to the contraction of the atomic 
Ad wavefunction from the left to the right of the periodic 
classification. 

In the d-dp model, the dd transitions are removed from 



Table VII: Hubbard-Kanamori parameters U,W,^7 (in eV) for 
the t2g-t2 g model and average interactions between the t2 s orbitals 
within d-dp in Sr2M04. The ratios U /V in t2g-t2g and U m m/vmm 
in d-dp have been multiplied by a factor 100. 





t2g-t2g 








d-dp 






(eV) 


u 


W 


J 


U/V 




Jm 


U mm J V mrn . 


Mo 


2.77 


2.15 


0.28 


24.8 


3.8 


0.48 


25.3 


Tc 


2.70 


2.07 


0.28 


24.3 


3.9 


0.52 


24.5 


Ru 


2.56 


1.94 


0.26 


23.2 


4.0 


0.55 


23.9 


Rh 


1.76 


1.18 


0.23 


16.6 


4.4 


0.61 


24.8 



the screening in order to obtain the Hubbard interac- 
tion matrices (Tab. VI and Appendix D). As the bare 
interaction, the on-site Hubbard interaction gets larger 
with the 4c? electron number (Fig. 7). This is rational- 
ized by the fact that the screening only slightly increases 
from Sr2MoC>4 to S^RhCU as indicated by the small de- 
crease of the ratio U mm /v mm (Tab. VII), from 25.8/100 
to 24.8/100. The increase of the screening is thus not able 
to counterbalance the effects due to the stronger orbital 
localization. 

The Hund's exchange interaction (Fig. 7) also repro- 
duces the atomic trend, increasing with the atomic num- 
ber. 



Hubbard parameters within the t2 g -t2 g 
Hamiltonian 



Alternatively, one can construct the low-energy Hamil- 
tonian that only includes the (d xy ,d xz ,d yz ) degrees of 
freedom. The energy windows Wt 2 are given in Tab. V. 

Since the pd charge transfer energy decreases through- 
out the series, the tail on the oxygen atomic sites of 
the downfolded local orbitals gets larger from Sr 2 Mo04 
to S^RhC^. Consequently, the orbital localization de- 
creases with the Ad electron number, in contrast to the 
atomic d wavefunctions. The trends for the bare interac- 
tions thus deviate from the atomic ones, which were pre- 
viously reported for the d-dp model (Fig. 7 and Tab. VII). 

The orbitally-resolved interactions (see Appendix D) 
calculated within the ti g -ti g Hamiltonian are more 
anisotropic than their analogues within the d-dp Hamil- 
tonian. This is another signature of the less spherical 
character of the downfolded orbitals within ti g -ti g - The 
largest interactions are interestingly obtained on the d xy 
local orbital, for both the screened and bare cases. 

The on-site ti g interaction li and Hund's exchange 
J exhibit trends that are similar to V and J7t>are 
(Fig. 7). The decrease from Sr2Mo04 to S^RhCU is even 
more pronounced. This is due to the screening which 
strongly and non-linearly increases between Sr2Ru04 and 
Sr 2 Rh04. Indeed, the ratio U /V within the ti g -ti g model 
(Tab. VII) is divided by almost a factor two through the 
series, whereas the ratio U mrn /v mrn is slowly decreasing 
within the d-dp model. 
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The difference in the screening between the two low- 
energy Hamiltonians comes from the transitions between 
t 2g and e g Kohn-Sham eigenstates, which are removed 
from the total polarization in the d-dp model. These 
transitions - by causing a notable increase of the screen- 
ing between Sr2Ru04 and Sr2Rh04 - are responsible for 
the smaller effective interactions in the latter. Within 
DFT-LDA (Fig. 6), the e s 's come closer to the Fermi level 
in the late Sr2M04 perovskites. The e g bands are even 
partially filled in Sr 2 Rh04 leading to a metallic screen- 
ing that contributes to lower the effective interactions. 
This is an artefact of the undistorted crystal structure 
of Sr 2 Rh04 and is not the case for the realistic distorted 
structure 45 . 



VI. PERSPECTIVES : USING INTERACTIONS 
CALCULATED WITHIN CRPA IN MANY-BODY 
CALCULATIONS 

A. Dynamical screening effects 

Both, calculations using a t 2g Hamiltonian and a d- 
dp Hamiltonian do exist in the literature for one of 
the compounds considered in this work, namely SrVC>3. 
Indeed, this compound has been chosen as a bench- 
mark system for nearly every new implementation of 
LDA+DMFT 14 ' 46 ' 47 ' 65,70 . In Rcf. [14] for example, such 
a calculation within a t 2g model was performed with a 
Hubbard-Kanamori U = 4 eV, yielding a quasi-particle 
renormalization of Z = 0.6. A second calculation used an 
energy window of -8.10 eV to 1.90 eV, spanning both, the 
bands used in the t 2g model and the oxygen-p dominated 
bands located between -8 and -2 eV, thus corresponding 
to the d-dp model in the notation of the present work. 
For the latter, a value of the intra-orbital interactions, 
Umm — 6.0 eV was used, leading to a similar quasi- 
particle renormalization (Z = 0.57) as in the t 2g -only 
model. 

The need for a larger interaction value in the d-dp 
model than in the t 2g model for reproducing the same 
mass renormalization is expected, and consistent with 
our results. The interaction values used in the calcula- 
tions of Rcf. [14] are however larger than what is found 
within cRPA: U = 3.2 eV for the t 2g model (Tab. Ill), 
Umm — 4 eV for the d-dp one (Tab. II). This fact cor- 
responds to the general observation of interaction val- 
ues used in many-body calculations in the literature be- 
ing quite often somewhat larger than the ones calculated 
from cRPA. 

This fact has remained a puzzle in the field for quite a 
while. Very recently a solution was proposed in Ref. [57], 
based on considerations including effects from the dy- 
namical screening of the Coulomb interactions in cor- 
related materials 55 ' 56 . These works concluded that the 
inclusion of dynamical screening effects resolves the ap- 
parent mismatch between cRPA values for the Hub- 
bard interactions and what is needed e.g. in standard 



LDA+DMFT to obtain agreement with experiments (e.g. 
to reproduce the experimental mass enhancement). We 
call here "standard LDA+DMFT" the usual procedure 
of supplementing a one-particle Hamiltonian obtained 
from LDA with static Hubbard (and Hund) interactions, 
and solving the resulting many-body Hamiltonian within 
DMFT. The extension of this procedure investigated 
in Refs. [55,56] takes into account also the frequency- 
dependence of the Hubbard U that results from the dy- 
namical nature of screening. To lowest order, the net ef- 
fect of this frequency-dependence is an additional renor- 
malization of the one-particle part of the Hamiltonian 
(see an explicit discussion of this effect in Ref. [57]). 
When this additional modification is taken into account 
- either implicitly by an explicit inclusion of dynami- 
cal U in the DMFT calculation, as done for SrVC>3 e.g. 
in Ref. [14], or explicitly by determining the additional 
renormalization factor from the frequency-dependence of 
U as in Ref. [57] - the correct mass renormalization is 
obtained, even when using the seemingly small ab ini- 
tio value obtained within cRPA. One can thus conclude 
that standard LDA+DMFT calculations need artificially 
increased U values to compensate for the lack of dynami- 
cal screening effects. If these effects are however included 
(e.g. by an explicit renormalization of the one-particle 
part of the Hamiltonian following Ref. [57]) the values 
calculated within the present work should be used. 



B. ti g -ti g or d-dp Models ? 

We have calculated, using our cRPA scheme, both, 
Hubbard and Hund interaction parameters suitable for 
t 2s -only and d-dp models. For materials where both 
choices yield a valid effective low-energy description of 
the physical properties, observables calculated within a 
many-body calculation using either model should give 
the same results. This has been verified explicitly for 
the case of SrVC-3 (e.g. in Ref. [14]), albeit with the 
caveat about additional renormalizations stemming from 
dynamical screening. Analogous results can be expected 
for early transition metal oxides, where hybridization ef- 
fects between d and p states are not too strong. This 
includes all the materials considered in the present work, 
except possibly SrMnC>3, where the high-spin character 
of the half-filled t 2g shell and the close-lying e g states 
may make a d-dp model more suitable for the calculation 
of some observables. 

In the general case, one has to be aware of a trade-off 
that has to be done in the construction of the low-energy 
model. Indeed, d-p hybridization is dealt with in a very 
different manner in the two model constructions that we 
have discussed. 

In the t 2g model, the Wannier functions arc built from 
the threefold degenerate low-energy bands only. While 
having majority t 2g character, these bands can contain 
substantial admixing of p-states, through the p-d hy- 
bridization. Downfolding then results in an "effective" 
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tig orbital, which becomes more and more extended the 
stronger the hybridization, displaying pronounced tails 
on the ligand sites. This may not be a problem as such, 
but the then quite extended "effective ti g " orbitals may 
make subsequent local approximations (in the construc- 
tion of the model, when inter-site interaction terms are 
neglected, or in the solution of the model, e.g. within 
DMFT) questionable. If the orbitals are too extended, 
one expects on the contrary also intersite interactions to 
become substantial, and the limitation to a Hubbard- 
type model with only on-site terms becomes question- 
able. 

At the same time, in late transition metal oxides, one 
encounters however a second, related problem that inval- 
idates the tig only model: the energetic position of the 
ligand p-states moves towards the Fermi level, and even- 
tually the gap forms not between Hubbard bands of t 2g 
character, but between ligand states and the upper Hub- 
bard band (charge transfer insulator) 64 . We have limited 
ourselves in this work to early transition metal oxides in 
order to investigate only cases where the i 2ff -only model 
makes physical sense. 

In the d-dp model, Wannier functions are built from 
both, bands with majority-<i and bands with majority-p 
character: the resulting Wannier functions are more lo- 
calized and thus more truly of ci-character than the ones 
in the t 2g model. Applying the local Hubbard interaction 
to the d states only is therefore a priori a well-defined pro- 
cedure. The question that arises in this context is how- 
ever the one about neglecting interorbital interactions, 
e.g. U p d between d and p orbitals, or U pp between the 
p-states. While neglecting the latter is likely a good ap- 
proximation because in reality the p-states are rather ex- 
tended and nearly full, the neglect of the former seems 
questionable when the p-d hybridization becomes strong. 
This is the reason for considering only early transition 
metal oxides in this work. We note that these questions 
are more questions about the nature of the appropriate 
low-energy model for a given compound than about cal- 
culating the parameters of such a model. Indeed, cRPA 
is well able to also yield U p d or non-local interactions, 
once the orbitals and subspaces are defined. 



C. Four-index U vs. Slater parametrization 

Following Eq. 16, a four-index U is calculated and 
could in principle be used in a many-body calculation 
without any restriction from the side of cRPA. Two-index 
or average interaction quantities are introduced only in 
the perspective of a further many-body calculation such 
as LDA+U or LDA+DMFT. In LDA+DMFT for exam- 
ple, even two-index spin-flip and pair-hopping terms (of 
energy scale J) are not routinely included due to the nu- 
merical cost of such calculation 43 . We note, for the case 
of SrV0 3 , that the largest three-index term is about 1CP 3 
eV, which compares to the smallest two-index one of 2.01 
eV (see Appendix C). This order of magnitude suggests 



that the more than two-index terms can reasonably be 
neglected. 

Of course, restricting the interactions to density- 
density (two-index) interactions does not necessarily im- 
ply that these have to be averaged. At this point, a subtle 
issue related to the construction of the one-particle part 
of the Hamiltonian comes into play. Indeed, LDA+U as 
well as LDA+DMFT comprise a "double counting cor- 
rection" meant to substract from the LDA Hamiltonian 
the mean-field like contribution due to interactions be- 
tween the correlated orbitals. This term amounts, in 
the standard formulation, to a global shift of correlated 
against uncorrelated orbitals. If the many-body interac- 
tions were very different for different correlated orbitals, 
one would expect that also the double counting would 
have to bear orbital-dcpcndcncc. This is an important 
but largely unexplored issue, and the standard proce- 
dure consists therefore in choosing spherically averaged 
(atomic-like) interactions. 

In practice, this question becomes of relevance only if 
the shape of the correlated orbitals is very anisotropic, 
so that the Slater parametrization becomes a poor ap- 
proximation. The quality of the latter is related to the 
degree of localization of the orbitals, so that the quality 
of the parametrization can to some extent be controlled 
by the choice of the energy window (the larger the energy 
window, the purer the atomic-like d-charactcr of the or- 
bitals). As we show in Appendix C, this parametrization 
is excellent for simple oxide materials such as SrVO,3. 



VII. CONCLUSIONS AND OUTLOOK 

Viewing the Hubbard and Hund interactions of a corre- 
lated electron material as bare interactions within a refer- 
ence system with reduced number of degrees of freedom, 
gives these parameters the status of auxiliary quantities, 
only meant to represent the physical fully screened inter- 
action W. From a conceptual point of view, they are then 
defined as soon as the reference system is specified. From 
a practical point of view, also the approximations to W 
and to the polarization of the reference system P sub have 
to be chosen. The lattice version of cRPA corresponds 
in this language to the choice of a low-energy subsystem, 
and the choice of the random phase approximation for 
the evaluation of all polarizations. 

In this work, we have presented an implementation of 
the constrained Random Phase Approximation in a den- 
sity functional theory electronic structure code within 
the linearized augmented plane wave framework. The 
method gives access to the matrix elements of the Hub- 
bard interaction matrix in a localized basis set of a 
downfolded lattice Hamiltonian. The strength of the 
Coulomb interactions is parametrized by the Hubbard 
U and Hund's exchange J, which can be used for the de- 
scription of correlated electron systems within interacting 
lattice Hamiltonians. 

We have calculated the Hubbard U and Hund's ex- 
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change J for the 3d and Ad ternary oxides SrMC>3 (M= 
V, Cr, Mn) and (M= Nb, Mo, Tc). Within a low-energy 
Hamiltonian ( "d-dp" ) that includes both d and oxygen p 
degrees of freedom - and an occupation energy cost on d 
only - U does not change much with the d electron num- 
ber. We have rationalized this trend as a competition 
between two aspects: increasing orbital localization - as 
evidenced also by the increasing bare interaction v - and 
increasing screening due to the reduction of the charge 
transfer energy with the oxygen ligands - as illustrated 
by a decrease of the screened interaction W. An alter- 
native low-energy Hamiltonian, based on the t2 g degrees 
of freedom only ("t2 g -t2 S ")> has also been constructed. 
Within this framework, both U and - to a lesser extent 
- J decrease with the d electron number in 3d oxides, as 
a consequence of the extended character of the localized 
t2 g orbitals. This is due to the downfolding of the tails 
of oxygen-p character. In Ad oxides, the trends for U and 
J are flat again, because of the larger charge transfer en- 
ergy with the oxygen-p states. The surprisingly smaller 
U for SrMn0 3 than for SrTc0 3 within cRPA for such 
Hamiltonian can be understood in terms of the weaker 
screening effects in the latter. 

We have considered analogous low-energy Hamiltoni- 
ans for the materials with the layered perovskite struc- 
ture Sr 2 M0 4 (M= Mo, Tc, Ru, Rh). For the same rea- 
sons as in the ternary oxides, U and J increase with the 
atomic number within the d-dp Hamiltonian, while they 
decrease within the £29-^29 Hamiltonian. 

We have emphasized the dependence of the effective 
Coulomb interactions on the choice of the one-particle 
part of the Hamiltonian. Determining Hubbard U and 
Hund's J entirely from first principles for a given low- 
energy Hamiltonian opens the way to a truly ab initio 
description of corrrelated materials within many-body 
calculations. 



ACKNOWLEDGMENTS 

We acknowledge useful discussions with M. Aichhorn, 
F. Aryasetiawan, M. Casula, A. Georges, D. Khom- 
skii, M. Imada, C. Martins, T. Miyake, J. Mravljc, L. 
Pourovskii, P. Rinke, G. Sawatzky, M. Scheffler, M. 
van Schilfgaarde and V. Vildosola. This work was sup- 
ported by the French ANR under project SURMOTT 
and IDRIS/GENCI under project 1393, and Natural Sci- 
ence Foundation of China (Project No. 20973009). 

While this work was being written, an alternative im- 
plementation of cRPA within a pseudopotential code was 
proposed in Ref. [103]. 



APPENDIX A : Technicalities of the cRPA 
implementation within the (L)APW+lo framework 

We report below on the technicalities that are specific 
to the cRPA part. Our cRPA implementation relies on 



the recent GW implementation within the (L)APW+lo 
framework called FHI-GAP 104 . Details related to the con- 
struction of an optimized product mixed basis for ex- 
panding polarizations, dielectric functions and Coulomb 
interactions, can be found in Ref. [104]. We will hence 
employ the same notations as in this reference. The prod- 
uct mixed basis, for example, will be denoted {x?( r )}> 
where (/-vectors belong to the first Brillouin zone and i- 
indices run over the size of the complete product mixed 
basis. 

The total polarization at the RPA level can be ex- 
panded into the product mixed basis {xf (?-)} 105 ' 106: 

P(r,r>;u) = ^EWWr p «(9.<( r ')> 

1 ij 

where the q-summation (of Af size) is performed over the 
first Brillouin zone. The polarization matrix elements in 
the product basis, Py, read as : 

occ.unocc 

Piifau,) = ^E E MU>(k,q)F nn , k (q,w)[M 3 nn ,(k,q)]*, 

k n.n' 

where the overlap between the Bloch wavefunctions, 
ipkn, and the product mixed basis functions, denoted 
M^ n ,(k,q), is defined as follows : 

M*„,(fc, 9 ) = / dr ^ kn (r)[ X f(r )1> k - q , n . (r)]*, 
and : 

^nk.n'k — q ^n'k — q ^nk 



U - 00 n k,n'k-q + 11] U + UJ n k,n'k-q - 11] 

Within the assumption that the correlated subspace 
C is unambiguously defined since the target correlated 
bands do not energetically overlap with the itinerant 
ones, it follows for the (i-restricted polarization P d (Eq. 9) 
expanded into the product mixed basis : 

occ,unocc 

Pt(q,u,) = ^E E M^(k,q)F dd , k (q,aj)[M^,(k,q)r, 

k d,d' 

and hence for the constrained polarization, P r : 
PTj(q,u) - Pij(q,u>) - I#(q,u,). 

The total symmetrized dielectric function, e, in the 
product mixed basis is defined as follows 104 : 

i i 
Sij (q, u) = % - v? (q)Pij (q, uj)v* (q) 

^ occ.unocc ^ 

= 6 V-jfH S 4{q)M i nn ,{k,q)F nn , k {q,oj) 

k n,n' 

x[vJ(q)Mi n ,(k,q)}*. 
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Analogously, the c?-constrained (partial) dielectric 
function, e r , reads as 



Sij 



(q)P[ j (q,Uj)vJ(q). 



From the Dyson-like equation for the screened inter- 
action (Eq. 5), it follows for W,W r expanded into the 
product mixed basis : 



part is expressed in terms of Slater integrals {F fc }. These 
Slater integrals are deduced from the interaction ma- 
trix elements computed by cRPA in the Wannier basis, 
{<^m-2<m<2} that is said "spheric" because of its com- 
plex representation 53 : 



F fc M = Ci, 



^ E 

mi, 7712, m.3, 7714 



C_1 \mi+m 4 7-r(spheric) / \ 



Wij(q,uj) 



v l 2 (q)e lJ 1 (q,uj)v J 2 (q) 



I k I 
-mi mi — 7713 m 3 



I k I 
-m 2 m 2 — mn 7714 



W^uj) = vf{qW{q,u)]^v]{q). 

The product Kohn-Sham representation of W r 
(Eq. 16) can be expanded into the product mixed basis 
via the overlap quantities M* n ,(fe, q) : 

(lpk ini 1pk2n2\W r \lpk 3 n 3 1pk 4 n 4 ) = ^ E [ M ™l«3 ( fel ' "? 

xW7 i (qr,w)M4„ 2 (fe 4 ,q)(5fc3,fc 1 _q(5 fe2!fe4 -q. 
and therefore, it follows for the Hubbard interactions : 



(35) 



where the parentheses correspond to the Wigner 3j- 
symbols and the coefficients C^k are defined as follows : 



2k + 1 



(21 + 1) 2 




(36) 



U 



Ft 1 FL2 R3 



YE"'" E^ n M 

1 ij 



xWT.( 9)W )M£ 4Jl4iLaHa ( 9 ), 
where M l LR LR , (q) are defined as follows : 



The frequency dependence of the Slater integrals arises 
from the frequency dependence of the interaction matrix 
elements induced by the dynamical screening. 



APPENDIX C : Reduced interaction matrices for 
SrMOu (M = V, Cr, Mn) 



M 



LR.L' 



We give below the reduced interaction matrices 
(Eqs. 19, 20 and 21) within the d-dp model (see Tab. I for 
the choice of the energy windows used to construct the 
d Wannier orbitals) and calculated with cubic symmetry. 
[P Ln (k, )]*M£ , (fe, q)PL' n ' (k — q)- In the following, the ordering of the orbitals in these ma- 
trices is d z 2 , d x 2 _ y 2 , d xy , d xz , d yz . The values are given in 
eV. 

In the case of SrVC>3, we also add the reduced inter- 
action matrices, U mm (Eqs. 22, 27, 28 and 29) that are 
deduced from the Slater integrals given in Tab. II. This 
allows for an estimation of the accuracy of the Slater 
parametrization. 



E 

n.n'EVi 



APPENDIX B : Parametrization of the Hubbard 
interaction matrix with Slater integrals 



The interaction matrix in atoms can be efficiently 
parametrized by a finite number of Slater integrals 23,24 ' 59 
thanks to the atomic sphericity. Such parametrization 
can be extended to solids if one assumes that within the 
solid : 



SrVOs 
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rr(sphcric) 
^ m\m2m 3 m 4 



^2a k (mi,m 27 m 3 ,m 4 )F k , (33) 



k=o 



where the angular part is described with Racah- Wigner 
coefficients. a^mi, m 2 , m 3 , 7714) are calculated with 
spherical harmonics, Yi m , as follows : 



a k (mi,m 2 ,m 3 ,m 4: ) 



q— — k 

x(Yi m2 Y kq \Yi m4 ), 



(34) 



((lj imi \Yi 2m2 Yi 3Tn3 ) corresponds to the Gaunt coefficient 
calculated with spherical harmonics), whereas the radial 



U° 



TTO-O- _ 





/ 4.43 


2.88 


2.73 


3.19 


3.19 \ 




2.88 


4.43 


3.35 


2.88 


2.88 




2.73 


3.35 


3.97 


2.75 


2.75 




3.19 


2.88 


2.75 


3.97 


2.75 




V 3.19 


2.88 


2.75 


2.75 


3.97 / 




f 


2.10 


2.01 


2.70 


2.70 \ 




2.10 





2.94 


2.24 


2.24 




2.01 


2.94 





2.15 


2.15 




2.70 


2.24 


2.15 





2.15 




^ 2.70 


2.24 


2.15 


2.15 


/ 
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Employing the Slater integrals F° = 3.2 eV, F 2 = 6.6 
and F 4 = 5.3 eV, we obtain for the Slater-parametrized 
reduced interaction matrices with cubic symmetry : 



U° 



Tract 

^ mm' 



u° 





/ a 1 n 
' 4.1U 


O 7Q 
Z. 16 


O 7Q 
Z. 16 


o.lo 


Q 1 Q \ 




Z.tS 


a 1 n 
4.1U 


Q QQ 
O.OO 


z.oo 


Z.OO 




Z. 1 6 


Q QQ 
O.OO 


A 1 ft 
4. 1U 


O QQ 


O QQ 
Z.OO 




o in 

o.lo 


O QQ 

z.ss 


O QQ 
Z.OO 


4.1U 


O QQ 
Z.OO 




V 3.18 


2.88 


2.88 


2.88 


4.10 / 




( ° 


2.01 


2.01 


2.68 


2.68 \ 




2.01 





2.91 


2.24 


2.24 




2.01 


2.91 





2.24 


2.24 




2.68 


2.24 


2.24 





2.24 




V 2.68 


2.24 


2.24 


2.24 


o / 




SrCrOz 








/ 3.84 


2.39 


2.40 


O Q7 
Z.O I 


O Q7 \ 
Z.O ( ' 




2.39 


3.84 


3.02 


Z.OO 


Z.OO 




2.40 


3.02 


3.89 


z.oy 


z.oy 




2.87 


2.56 


2.59 


Q QO 

o.oy 


z.oy 




V 2.87 


2.56 


2.59 


2.59 


3.89 / 




( ° 


1.66 


1.67 


2.38 


2.38 \ 




1.66 





2.61 


1.91 


1.91 




1.67 


2.61 





1.94 


1.94 




2.38 


1.91 


1.94 





1.94 




V 2.38 


1.91 


1.94 


1.94 


/ 




SrMnOs 








/ 3.62 


2.17 


2.25 


2.74 


2.74 \ 




2.17 


3.62 


2.90 


2.41 


2.41 




2.25 


2.90 


3.90 


2.52 


2.51 




2.74 


2.41 


2.52 


3.91 


2.52 




V 2.74 


2.41 


2.51 


2.52 


3.90 / 



APPENDIX D : Reduced interaction matrices for 
Sr 2 M0 4 (M = Mo, Tc, Ru, Rh) 

u d-dp Hamiltonian" 

The reduced interaction matrices within the d-dp 
model (Tab. V) and calculated with cubic symmetry, are 
given below for the layered perovskites. The ordering of 
the orbitals in these matrices is d z 2,d x 2_ y 2,d xy , d xz ,d yz . 
The values are given in eV. 



Sr 2 MoOi 



U" 





( 4.20 3.04 


2.90 


3.22 


3.22 \ 




3.04 4.35 


3.43 


3.06 


3.05 




2.90 3.43 


3.97 


2.93 


2.92 




3.22 3.06 


2.93 


3.86 


2.89 




V 3.22 3.05 


2.92 


2.89 


3.84 J 




/ 2.48 


2.38 


2.83 


2.83 \ 




2.48 


3.10 


2.56 


2.55 




2.38 3.10 





2.44 


2.44 




2.83 2.56 


2.44 





2.41 




\ 2.83 2.55 


2.44 


2.41 


/ 


Sr 2 Tc0 4 








/ 4.06 2.88 


2.81 


3.13 


3.13 \ 




2.88 4.23 


3.38 


2.96 


2.96 




2.81 3.38 


4.04 


2.90 


2.90 




3.13 2.96 


2.90 


3.86 


2.84 




V 3.13 2.96 


2.90 


2.84 


3.86 / 




/ 2.29 


2.25 


2.72 


2.72 \ 




2.29 


3.04 


2.44 


2.44 




2.25 3.04 





2.38 


2.38 




2.72 2.44 


2.38 





2.33 




V 2.72 2.44 


2.38 


2.33 


/ 



Sr 2 RuOi 



( 


1.44 


1.49 


2.22 


2.22 \ 


1.44 





2.47 


1.73 


1.73 


1.49 


2.47 





1.81 


1.81 


2.22 


1.73 


1.81 





1.81 


V 2.22 


1.73 


1.81 


1.81 


J 



XJ aa , = 



( 4.06 2.83 2.78 3.17 3.18 \ 

2.83 4.22 3.37 2.98 2.99 

2.78 3.37 4.07 2.94 2.95 

3.17 2.98 2.94 4.02 2.93 

V 3.18 2.99 2.95 2.93 4.05 / 
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u: 



era 
mm' 



^ mm' 



( 2.22 2.19 2.73 2.74 \ 

2.22 3.01 2.42 2.44 

2.19 3.01 2.39 2.39 

2.73 2.42 2.39 2.38 

\ 2.74 2.44 2.39 2.38 

Sr 2 Rh0 4 



U° 



U acr , = 



( 4.18 2.93 2.95 3.37 3.38 \ 

2.93 4.34 3.55 3.15 3.16 

2.95 3.55 4.37 3.18 3.19 

3.37 3.15 3.18 4.39 3.20 

\ 3.38 3.16 3.19 3.20 4.44 / 



Sr 2 TcOi 



( 2.89 2.10 2.10 \ 

2.10 2.61 2.01 

2.10 2.01 2.61 

0.00 1.81 1.81 

1.81 0.00 1.74 

\ 1.81 1.74 0.00 J 



Sr 2 Ru0 4 



U° 



( 2.29 2.31 2.92 2.94 \ 

2.29 3.17 2.56 2.58 

2.31 3.17 2.57 2.58 

2.92 2.56 2.57 2.59 

\ 2.94 2.58 2.58 2.59 / 



2.72 1.97 1.97 
1.97 2.48 1.89 
1.97 1.89 2.48 



il t2g-t2 g Hamiltonian" 



We give below the reduced interaction matrices for 
the tig local orbitals within the t 2g -t2 g model. The 
Hubbard-Kanamori parameters, U and U\ shown in Ta- 
ble VII correspond to the orbital average, U™ m ,XJ%L. m , 
(Eqs. 19 and 21), respectively. The exchange param- 
eter, J , corresponds to the orbital average of J m m' — 
(U™ m ,-U™ m ,)(l-5 mm ,) (Eq. 20). The ordering of the 
orbitals is (d xy , d xz , d yz ). 



I 0.00 1.71 1.71 \ 

1.71 0.00 1.65 
v 1.71 1.65 0.00 j 



Sr 2 Rh0 4 



Sr 2 MoOi 



2.96 2.19 2.19 

U™ m , = [ 2.19 2.68 2.09 

2.19 2.09 2.68 

0.00 1.89 1.89 

/ = I 1.89 0.00 1.82 

1.89 1.82 0.00 



u° 



( 1.81 1.16 1.20 
1.16 1.69 1.19 
v 1.20 1.19 1.77 



/ 0.00 0.93 0.96 \ 

0.93 0.00 0.97 
v 0.96 0.97 0.00 j 
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